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^ . Abstract. 

O i BeppoSAX observed the z=3.27 quasar PKS2126-158 on 1999 May 24-28 when its 2-lOkeV and 0.1-2.5keV 

CO . fluxes were 1.1x10"" and 4.4x10"^^ erg cm s respectively, a factor of 2 higher than in all previous ROSAT 

'/^ ' and ASCA observations and 40% higher than in two more recent Chandra and XMM-Newton observations. The 

shortest detected rest frame variability timescale is of a few months, comparable to the causal timescale associated 

■ to an emission region of ~ 10 Schwarzschild radii around a fewxlO^^M© black hole. The source is detected with 
a signal to noise ratio up to ~50 keV, 215 keV rest frame. The BeppoSAX observations confirm the 
presence of low energy absorption along the line of sight, independent on the continuum model adopted, at high 

• confidence level. Despite the limited spectral resolution of the BeppoSAX LECS and MECS it is possible to put 

constraints on different absorption and continuum models, but not to unambiguously determine the redshift of 
the absorber. If the absorber is not significantly ionized the BeppoSAX data do prefer an absorber at z < 2.7. 
5_j ' Strong and complex metal line systems along the line of sight to PKS2126-158 have been found at z=0.6631 and 

^ I at 2.64<z<2.82. They could well be associated to the X-ray absorption. Conversely, an ionized ("warm") absorber 

. at the quasar redshift provides a good fit only if the iron abundance is smaller than ~ 0.3 solar, while that of 

the other elements is fixed to the solar value. Lower iron metallicity would imply a lower dust to gas ratio, since 
iron aggregates easily in dust. This can help in solving the apparent paradox of the lack of significant ultraviolet 

■ reddening in this source while strong absorption is detected in X-rays. Low iron abundance would be at odds with 
$H ' the supersolar abundances derived from the broad emission lines. 

Key words, quasars spectrum - quasars high redshift - quasars - metal abundances 



1. Introduction 

The high redshift (z>3) Universe can be probed by observing distant convenient "lighthouses" such as quasars, Gamma 
Ray Bursts, and luminous galaxies. Both emission and absorption features superimposed onto the source continuum 
can be used for this purpose. High quality spectra of high redshift cosmic sources can tell us about the physics and 
chemistry of the matter along the line of sight (i.e. density, physical state and element abundances). In particular, 
high-z quasars and galaxies give a unique chance of probing the metallicity of the gas enriched by the first generations 
of stars. Iron is expected to be produced mostly by type la SNe, occurring with a delay of « 1 Gyr with respect to 
the first type II SNe, which produce most of the a elements such as O, Si, Mg. Therefore, the relative abundances of 
various elements produced on different timescale in high-z quasars (e.g. the [a/Fe] logarithmic ratio, where the columns 
are normalized to their solar value) is a sensitive clock of star-formation history, and its determination provides unique 
constraints on the star-formation and on the initial mass function at high redshift (e.g. Hamann & Ferland 1999). 

So far metallicities at high-z have been studied by means of a few quasar and radio-galaxy UV rest-frame emission 
lines (see e.g. Hamann & Ferland 1999, Constantin et al. 2002, Hamann et al. 2002), and via optical-UV absorption 
systems, both associated with the quasar and with galaxies which happen to lie along the line of sight to bright quasars. 
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Emission line metallicity determinations are plagued by uncertainties on the properties of the emitting gas (tem- 
perature, density, ion fraction, dust depletion) and by the small number of lines available to constrain the metallicity. 
In particular, the abundance of iron relative to the a elements is mostly inferred from the intensity of the UV Fe 
bump relative to Mgll (2803A) (Thompson et al. 1999, Iwamuro et al. 2002, Dietrich et al. 2002). However, the Fell 
UV bump is one of the most important coolants of the transition region of the BLR in AGNs and therefore might be 
characterized by a sclf-rcgulating intensity; indeed, Verner et al. (1999) have shown that the emissivity of the Fc UV 
bump is far from having a 1:1 relation with the iron abundance. As a consequence, the iron metallicity of high-z QSOs 
inferred by a naive comparison of the Fell-UV/Mgll intensity at high and low redshift might be quite inaccurate or 
even deceiving. 

Quasar intrinsic broad absorption line (BAL) metallicity determinations provide bizarre but usually supersolar 
metal abundances (Hamann 1998, Hamann & Ferland 1999, Arav et al. 2001, Hamann et al. 2003). However these 
could have be overestimated, because strong absorption lines like CIV could be optically thick and strongly saturated 
(e.g. Hamann et al. 2003). The analysis of BALs is complicated by the fact the usually the line troughts do not reach 
zero intensity, because the absorbers are leacky and/or because of scattered/emitted flux from the absorbing region. 
Furtheremorc, from the optical/UV BAL it is not possible to measure the abuncances of iron, which is one of the key 
elements to constrain the star formation history. 

In the case of absorption features not associated with the quasars, optical-UV metal line systems in galaxies along 
the line of sight to quasars probe mainly the ISM of outer haloes, rather than the bulge or the disk (see e.g. Churchill 
et al. 2000). Even the highest column density systems known, the Damped Lya systems (DLAs), may give a distorted 
view of the bulk of the 'typical galaxy' population at z=2-4. DLA studies are likely to be biased against dusty systems, 
since these would hide the background quasars in color-based surveys (e.g. Savaglio et al. 2000). Furthermore, in 
DLAs metallicities are usually evaluated using faint zinc lines, since iron, which is 20 times more abundant than zinc, 
is heavily depleted onto dust grains (see e.g. Pettini et al. 1997). Using zinc lines Pettini et al (1999) and Prochaska 
& Wolfe (1999) estimate a metallicity of about 0.1 at z=3-4. DLA metallicities may approach solar values at low-z 
(Savaglio et al. 2000). 

Summarizing, metallicity measurements at high redshift based on optical observations have often provided only 
partial or uncertain information. This statement applies specifically to iron. In contrast. X-ray measurements offer 
a unique means of investigating absorption systems because the X-ray absorption edges can constrain tightly the 
ionization state, the molecular state and the dust content of intervening material and so probe total column densities. 
A direct measurement of the depth of iron X-ray absorption edges will remove the degeneracy involving various physical 
parameters of the emitting region, which plague current metallicity estimates. This can be achieved at z=3-4 with 
the present instrumentation, because Fe-K features at these redshift are redshifted to 1.5-2 keV, where the X-ray 
instrument sensitivity is the greatest. 

This approach is plausible as there are strong signs of absorption toward several z=2-3 radio-loud quasars with 
implied column densities of Nh ~ 10^^ cm~^ in the quasar rest frame (Elvis et al. 1994, Serlemitsos et al. 1994, Cappi 
et al. 1997, Fiore et al. 1998, Reeves & Turner 2000). The quality of the ROSAT and ASCA discovery spectra however 
can not discriminate between a location for the absorber at the quasar or along the intervening line of sight. Fiore 
et al. (1998) and Elvis et al. (1998) argue for intrinsic origin because of the lack of absorbers in radio-quiet quasars 
at z=2-3. Kuhn et al. (2001) show that z=3 radio-quiet quasars have spectra that are as UV-bright as low redshift 
radio-quiet quasars, implying Ay ^ 0.3 mags. However the absence of dust and of detectable X-ray Nh in z=2-3 radio 
quiet quasars could be a selection effect. Most quasars are selected via optical colors and a fairly modest reddening can 
drop the objects out of the samples. Furthermore, radio-quiet quasars are fainter in X-rays than radio-loud quasars 
and so have more poorly constrained spectra. So the issue of the location of the absorbers in radio-loud quasars is still 
open. 

PKS2126-158 is the brightest quasar in X-rays at z> 3, and it is the second brightest, after PKS2149-306, at z> 2. 
For this reason PKS2126-158 has been studied in great detail at X-ray frequencies since the first Einstein detection 
(Worrall & Wilkes 1990). Elvis et al. (1994) discovered a strong low energy cut-off, which is naturally explained by 
photoelectric absorption along the line of sight, in addition to the Galactic 21cm column density of 4.85 x 10^*^ cm~^ 
(Elvis et al. 1989). The implied column density ranges from 0.9 — 20 x 10^^ cm^^, if the absorber is at 0, it is not 
highly ionized and has solar abundances, up to a column of 0.8 — 2.7 x 10^^ cm~^ if the absorber is at z~ 3.27, the 
quasar redshift. ASCA observations confirmed the low energy cutoff (Serlemitsos et al 1994, Cappi et al. 1997, Reeves & 
Turner 2000) but cannot distinguish beween genuine photoelectric absorption or continuum curvature. We have taken 
advantage of the good BeppoSAX sensitivity over a very broad band, from 0.1 keV up to 200 keV, to disentangle 
photoelectric absorption and continuum curvature and to constrain the redshift, ionization state and metallicity of the 
absorber in PKS2126-158. 

The paper is organized as follows: Section 2 presents the observations and data reduction; Section 3 presents the 
spectral analysis; Section 4 discusses the main results on the low energy absorption. 
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2. BeppoSAX data reduction and analysis 

BeppoSAX observed PKS2126-158 from 1999 May 24 through 1999 May 28. The observations were performed with the 
BeppoSAX Narrow Field Instruments, LEGS (0.1-10 keV, Parmar et al. 1997), MEGS (1.3-10 keV, Boella et al. 1997b), 
HPGSPC (4-60 kcV, Manzo ct al. 1997) and PDS (13-200 kcV, Frontcra ct al. 1997). LEGS and MEGS arc imaging gas 
scintillation proportional counters, the HPGSPG is a coUimated high pressure gas scintillation proportional counter 
and the PDS consists of four phoswich units. The PDS is operated in the so called "rocking mode", with a pair of units 
pointing to the source while the other pair monitor the backgroimd ±210 arcmin away. The units on and off source 
are interchanged every 96 seconds. We report here the analysis of the LEGS, MEGS and PDS data; the HPGSPG data 
are very noisy, due to the high HPGSPC background. The MEGS observations were performed with units 2 and 3 
(on 1997 May 6*^ a technical failure caused unit MEGSl to be switched off); these data were combined together after 
gain equalization. The LEGS was operated during dark time only, therefore LEGS exposure times are usually smaller 
than MEGS ones. Table 1 gives the LEGS, MEGS and PDS exposure times and the count rates. 

Standard data reduction was performed using the SAXDAS software package version 2.0 following Fiore, Guainazzi 
and Grandi (1999). In particular, data are linearized and cleaned of Earth occultation periods (we accumulated data 
for Earth elevation angles > 5 degrees) and unwanted periods of high particle background (satellite passages through 
the South Atlantic Anomaly and periods with magnetic cut-off rigidity > 6 GcV/c). The low and stable particle 
background of the LEGS, MEGS and PDS (with variations at most of 30% around the orbit, due to the low inclination 
of the satellite orbit, 3.95 degrees), ensures that systematic errors in the subtraction of the background are negligible, 
i.e. smaller, in the energy bands of interest, than the statistical errors. 

LEGS and MEGS spectra were extracted from regions of 8 arcmin and 3 arcmin radii respectively. These radii 
maximize the signal-to-noise ratio below 1 keV in the LEGS and above 2 keV in the MEGS. Background spectra were 
extracted in detector coordinates from high Galactic latitude 'blank' fields (98_11 issue) using regions equal in size to 
the source extraction region. We have compared the mean level of the background in the LEGS and MEGS "blank 
fields" observations to the mean level of the background in the PKS2126-158 observations using source free regions at 
various positions in the detectors. The "local" MEGS background count rate is within 2% that in the "blank fields" . 
On the other hand, for the LEGS the "local" background count rate is 10% smaller than that of the "blank fields" . 
For this reason, the "blank fields" background spectrum was multiplied by a factor of 0.9 before subtraction from the 
source spectrum. 

The PDS data were reduced using the "variable risetime threshold" technique to reject particle background (see 
Fiore, Guainazzi & Grandi 1999). This technique reduces the total 13-200 keV background to about 20 counts s~^ 
and the 13-80 keV background to about 6 counts s~^ (instead of 30 and 10 counts s~^ respectively, obtained using the 
standard "fixed risetime threshold" technique). 

The PDS rocking mode provides a reliable background subtraction. This can be checked by looking at the spectrum 
between 200 and 300 keV, where the effective area of the PDS to X-ray photons is small and therefore the source 
contribution is negligible. After background subtraction we obtain a count rate of -0.0014±0.0100 counts s^^, fully 
consistent with the expected value of 0. The net (background subtracted) 13-100 keV on-source signal is 0.157±0.021 
counts s^^. The source is detected with a signal to noise ratio S/N^ 3 up to '--^50 keV, or 215 keV in the rest frame). 
Gonfusion in the PDS collimator Field Of View (1.4 degrees FWHM) ultimately limits our capability to constrain the 
high energy spectrum. Wc have carefully checked for any possible contaminant in a region of 1.5 degrees radius around 
the source (using the NED, SIMBAD, AGN, clusters, CVs. Radio and X-ray sources catalogs) finding no obvious 
bright hard X-ray source. Of course there is the possibility of a bright "unknown" source in the PDS field of view. 
The chance of finding a source in any given 2 square degrees, the PDS beam area, is however small. The HEAO-1 A4 
all sky catalog (Levine ct al. 1984) lists just 7 high Galactic latitude sources in the 13-80 keV band down to a flux of 
2x10^^" erg cm~^ s~^ (lOmGrab). The 13-80 keV flux is about 20 times smaller than this figure and so, assuming a 
logN-logS slope of —1.5, we expect a chance coincidence rate of 2%. 



Table 1. Observation log 



Instrument 


Exposure (ks) 


Count rate " 


LEGS 


78.1 


0.1-4 keV 5.50±0.09 


MEGS 


107.9 


1.7-10 keV 11.9±0.10 


PDS 


51.6 


13-100 keV 15.7±0.21 


" 10"^ counts 


s-i 





Spectral fits were performed using the XSPEG 11.0.1 software package and public response matrices as from the 
1999 December release. For the MEGS and PDS we used the standard on axis matrices (effective area and response 
matrix). For the LEGS we built the efl[ective area file corresponding to the source position in raw coordinates, since it 
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data and folded model 
PKS8136-15B BeppoSAX LEGS, MECS, PDS 




0.1 1 10 100 

channel energy (keV) 



Fig. 1. The LEGS, MECS and PDS spectra of PKS2126-158 fitted with a power law model with Galactic absorption. 

turned out that this position is relatively far from the default pointing position (by about 3 raw pixels, i.e. 45 arcsec.) 
This caused the source to be closer than usual to one of the tungsten wires supporting the LEGS window. Our LEGS 
effective area file takes than correctly into account the photons absorbed by the wire. 

LEGS and MEGS spectra were rebinned following two criteria: a) to sample the energy resolution of the detectors 
with four channels at all energies whenever possible, and b) to obtain at least 20 counts per energy channel. 

Gonstant factors have been introduced in the fitting models in order to take into account the intercalibration 
systematics between the instruments (Fiore, Guainazzi and Grandi 1999). Assuming the MEGS as reference instrument, 
the expected factor between LEGS and MEGS is about 0.9 [0.7-1.1]. Fitting the LEGS and MEGS spectra in the 
common 1.7-4 keV band with a simple model produced a factor of 0.835±0.025. In all following fits the LEGS-MEGS 
factor is therefore constrained to vary between 0.81 and 0.86. 

The expected factor between the PDS and MEGS is 0.8. We constrained the PDS normalization factor to vary in 
the range 0.7-0.9, to allow for confusion in the PDS on-source and off-source collimators due to faint hard sources. 
The energy range used for the fits are: 0.1-4 keV for the LEGS (channels 12-400), 1.65-10 keV for the MEGS (channels 
37-220) and 13-200 keV for the PDS. Errors quoted in this paper are 90% confidence intervals for two interesting 
parameters, unless differently specified. 

3. Results 

Figure ^ shows the LEGS-f MEGS-f PDS spectra fitted with a power law model absorbed by the Galactic column 
density along the Hue of sight (4.85 x 10^° cm~^ ^ Elvis et al. 1989). The figure shows that the 0.1-100 keV spectrum 
of PKS2126-158 is complex. A strong cut-off below 1 keV is clearly visible as well as a broad-band spectral curvature. 

To allow an easier comparison with results obtained with other satellites we give in table 2 the results of a fit 
with a simple power law model reduced at low energy by the Galactic column density of 4.85 x 10^" cm~^ and by an 
additional column at zero redshift, assuming solar abundances. The is acceptable, but relatively large residuals are 
present at both low (0.4-0.8 keV) and high (> 20 keV) energies. Both the spectral index and the absorbing column 
are consistent with previous ROSAT (Elvis et al. 1994) and ASGA (Serlemitsos et al. 1994, Gappi et al. 1997, Reeves 
& Turner 2000) determinations. 

Allowing the redshift of the absorber to be free to vary does not improve the fit. Figure [3 shows the confidence 
contours for the Nh and the redshift for a power law model continuum. With this choice for the continuum shape we 
can conclude that absorption in addition to the Galactic value is required by the data at high confidence level and 
that the redshift of the absorber must be smaller than 0.8 at the 99% confidence level. 

However, figure ^ indicates some curvature of the underlying continuum. This was already suggested by ASGA 
observations (Serlemitsos et al. 1994, Gappi et al. 1997). Intrinsic curvature of the low energy spectrum will of course 
result, in low resolution spectra, in a change of the shape of the actual low energy photoelectric cutoff and therefore 
of the estimated absorber properties (density, redshift, physical and chemical composition). 

We have therefore fitted the data using a power law with a high energy exponential cutoff (see table 2). The 
improvement in x^ with respect to the simple power law model is significant at > 99.99% confidence level. The ratio 
between the data and the best fit model is shown in figure El Allowing the redshift of the absorber be free to vary does 

^ Cappi et al. (1997) verified that large molecular clouds are not present along the line of sight to PKS2126-158, and that the 
Galactic gas column density estimated from the IRAS 100^ maps, assuming a typical gas to dust ratio, is ^ than the 21cm 
estimate. 
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Confidence conto 
PKS2136-158 BeppoSAX LECS, MECS, PDS 




nH 10^22 



Fig. 2. Nh vs redshift 67%, 90% and 99% confidence contours for a power law continuum. 



data/model 
PKS2126-158 BeppoSAX LECS, MECS, PDS 




channel energy (keV) 



Fig. 3. Ratio between the BeppoSAX LECS, MECS and PDS spectra of PKS2126-158 and the best fitting power law 
with exponential cutoff model. 

Table 2. Spectral fits witii a neutral absorber 



model 




Z 






^break Or Ec„t 


X'(dof) 


PL 


0.16±0.03 


OFIX 


0.62±0.03 






145.9/144 


CUTOFFPL 


Q ,0+0.05 


OFIX 


47+0-" 

-0.07 






131.8/143 


BKNPL 


lS+0 04 
^•r'J_o.03 


OFIX 


n 70+0.32 
''J-O.OS 


1.52±0.08 


5.oi?:° 


136.1/142 


Curved 


^•^■3-0.03 


OFIX 


n 4q+0 09 

"•^^-0.26 


r+3.0 
■-'■'^-1.8 


36+1? 


131.8/142 



" in 10^2 cni-2; 
power law energy index, F{E) oc E~"-'^; 

hard power law energy index in the BKNPL and Curved models, F{E) oc E^"'^ for E < E^reak and F{E) oc 



E-"" X eI'^^^i^""'' for E > Ebreak- The Curved model is similar to the BKNPL model but the two power laws are 
joined smoothly; 

Break energy in keV in the BKNPL and Curved models; 
" Exponential cut-off energy in keV in the CUTOFFPL model, F{E) oc £'~"^ x e^/^™'. 



not improve the fit. Figure^ shows the Nh vs. redshift confidence contours for the power law with exponential 
cutoff model continuum. 

We have also fitted the data using a broken power law model and with with a broken power law model with a 
smooth conjunction between the two power laws ("Curved" model). The results are again in Table 2. For the broken 
power law model The improvement in with respect to the simple power law model is significant at the 99.28% 
confidence level. Figure |31 shows the Nh vs. redshift confidence contours for a broken power law model continuum. 
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Confidence contours 
PKS3126-158 BeppoSAX LECS. MECS, PDS 




nH 



Fig. 4. Nh vs redshift 67%, 90% and 99% confidence contours for a broken power law continuum. 



Confidence contours 
PKS2126-158 BeppoSAX LECS. MECS. PDS 




Fig. 5. Nh vs redshift 67%, 90% and 99% confidence contours for a power law with exponential cutoff continuum. 



There are two minima in the surface, one at z—0 and the second at 2.2. The best fit results for the "curved" 
model resemble those of the power law with exponential cutoff model, but with an additional free parameter. 

Figures 01 and show that if the absorption is neutral an the continuum is paremeterized with either a broken 
power law or a power law with a cutoff, the redshift of the absorber is less that that of the quasar at better than the 
99% confidence level. 

For the three continuum models, we have also fitted an ionized absorber model (Absori in XSPEC) at both 
z=0 and z=3.27. The results are given in Table 3. The continuum parameters are very close to those in Table 2 and 
therefore, for sake of simplicity, we give in Table 3 the parameters of the absorber only. The fits with the absorber 
at z=0 give a good in cases. Equally good fits for an absorber at the quasar redshift can only be obtained if 
the iron abundance Zpe is substantially less than solar. (The iron K edges are redshifted into the 1.6-2.0 keV band. 
At these energies the quality of the LECS and MECS spectra can strongly constrain the iron abundance.) Since the 
abundance of all other elements is fixed to the solar value, the fits actually constrain the [a/Fe] ratio to be higher than 
0.5. 

As an example Figure shows the confidence contours for the iron abundance versus the ionization parameter, 
for a power law continuum model, assuming that the absorber is at the redshift of the quasar. The analogous contours 
for the other two parameterizations of the continuum are similar. In all cases the best fit solution implies column 
densities of the order of a few x 10^'^ cm"'^, high ionization parameter and low iron abundance. The inclusion of an 
additional cold absorber at the quasar redshift does not improve the x^- 
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Table 3. Spectral fits with an ionized absorber' 



model 




z 






x"(dof) 


PL 


0.26±0.05 


OFIX 




IFIX 


134.7/143 


BKNPL 


n 90+0.20 
"■^^-0.10 


OFIX 


< 0.8 


IFIX 


131.4/140 


CUTOFFPL 


20+°-^° 


OFIX 


< 1.0 


IFIX 


126.5/142 


PL 




3.27F1X 


18001-™ 


< 0.25 


140.8/142 


BKNPL 




3.27F1X 


1 300+=*™° 

iOUU_j]^00 


< 0.25 


129.5/141 


CUTOFFPL 


36+lf 


3.27F1X 


1700+=***°° 


< 0.3 


127.6/141 



^ continuum fit values are indistinguishable from those in Table 2; 

^ in 10^^ cm~^; ionization parameter; ^ iron abundance in unit of the solar value. 

Confidence contours 



^ PKS2126-158 BeppoSAX LEGS, MEGS, PDS 




1 0.2 0.3 0.4 

Fe_abund 



Fig. 6. Ionization parameter vs. iron abundance 67%, 90% and 99% confidence contours for a power law with 
exponential cut-off continuum. 



4. Discussion 

4.1. Flux, Luminosity and variability 

The 2-10 keV and 0.1-2.5 keV fluxes observed by BeppoSAX are 1.1 x 10^^^ and 4.4 x 10~^^ erg cm^^ s^^ respectively. 
The unabsorbed 0.1-2.5 keV flux is 2.1 times higher (1.0 x 10~^^ erg cm~^ s~^ ). The statistical uncertainties on 
these values is of the order of a few %, smaller than the systematic uncertainty in the LEGS and MEGS flux absolute 
calibration (which is of the order of 5 %, Sacco, private communication, also see Fiore, Guainazzi & Grandi 1999). 
The unabsorbed fluxes correspond to a 2-10 keV and 0.1-2.5 keV luminosities of 1.6 x 10** erg s~^ and 1.5 x 10** 
erg s~* respectively, assuming "isotropic" emission and using Ho=50 and ilo=0. For a Ho=65, rio=0.3 and f^A — 0.7 
cosmology the 2-10 keV and 0.1-2.5 keV luminosities are 2.5 times smaller (6.8 x 10*^ erg s^* and 6.2 x 10*^ erg s~* 
respectively) . 

The observed fluxes are compared in flgureQwith previous measurements obtained with the Einstein, EXOSAT, 
ROSAT and ASGA satellites (WorraU & Wilkes 1990, Elvis et al. 1994, Gappi et al. 1997) and with the fluxes recorded 
in two recent Ghandra and XMM-Newton observations (D'elia et al. 2003 in preparation, Ferrero & Brinkmann 2003). 
5% systematic errors have been added in quadrature to the statistical errors on the fluxes. Figure [3 shows that 
BeppoSAX caught the source in its historical maximum, more than a factor of 2 brighter than all previous observation 
and 40% brighter than in the later Ghandra and XMM-Newton observations. A detailed analysis of the source spectral 
variability among the ASGA, BeppoSAX, Ghandra and XMM-Newton observation is beyond the scope of this paper 
and will be presented elsewhere (D'Elia et al. 2003). The rest frame time delay between the ASGA and BeppoSAX 
observations is about one year, while that between the BeppoSAX and Ghandra and BeppoSAX and XMM-Newton 
observations is of 1.5 and 6 months respectively. 

PKS2126-158 is a GigaHertz Peaked Spectrum (GPS) quasar (O'Dea et al. 1991, Stanghelhni et al. 2001). On a 3 
miliarcsec scale the source has core-jet morphology, with a possible counter jet detection. Based on the possible counter 
jet detection and, more importantly, on the low polarization, Stanghelhni et al. (2001) conclude that PKS2126-158 can 
be only marginally beamed. If beaming is not large in this source, then the above "isotropic" luminosities are likely 
to be close (within a factor of a few) to the real luminosities. These then would imply extreme black hole masses of 
at least a fewxlO^^M©, assuming that the quasar is emitting close to its Eddington luminosity. The causal timescale 
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Fig. 7. 2-lOkeV (filled circles) and 0.1-2keV (open triangles) fluxes observed by Einstein, EXOSAT, ROSAT, ASCA, 
BeppoSAX and XMM-Newton in tlie last 22 years. Error bars are 1 a statistical uncertainties 



associated to an emission region of ~ 10 Schwarzschild radii around these black holes is ^months, comparable to the 
shortest variability timescale detected. 

The absorbing column density depends of course by the redshift. If z=0 is assumed then the BeppoSAX values in 
Table 2 are well consistent with the ROSAT and ASCA ones. If z—3.27 is assumed then the column densites in figures 
Eland 111(2 - 4 x 10^^ cm~^), are somewhat higher than the ROSAT and ASCA values (Elvis et al. 1994, Cappi et 
al. 1997), but still statistically compatible with them. We note that this little discrepacy is likely due to the fact that 
the continuum assumed by the latter authors is a simple power law, while figures ^ and are obtained using a broken 
power law and a power law with exponential cutoff continuum. 



4.2. Low Energy absorption 

The BeppoSAX NFI observations of PKS2126-158 confirm the presence of low energy absorption in this source indepen- 
dent of the continuum form. Unlike all previous investigations, even assuming a curved continuum and a broken power 
law continuum, the spectra do require, at a high confidence level, absorption in excess to the Galactic column density 
along the line of sight. Because of the limited spectral resolution of the LECS and MECS it is difficult to distinguish 
between different absorption and continuum models and to assess precisely the redshift of the absorber. If a power law 
model is used to parameterize continuum the redshift of the absorber is smaller than 0.8 at the 99% confidence level. 
If the continuum is parameterized using a broken power law the redshift of the absorber is more poorly constrained 
(the 99% contours nearly touch the emission redshift of 3.27). Interestingly there are two minima in the surface, 
one at z=0 and the second at z=2.2. If the continuum is parameterized using a power law with an exponential cutoff 
then Zabs < 3 at the 99% confidence level. In summary, if the absorber is not strongly ionized, the BeppoSAX data 
do prefer an absorber at z<Zem- X-ray absorption up to a level of 10^^ cm~^ can be produced by DLA systems along 
the line of sight. This is seen in several low-z DLAs by Chandra (Bechtold et al. 2001, Siemiginowska et al. 2002). 
Although, no DLA is observed toward PKS2126-158, only a small part of the redshift range is accessible to direct DLA 
searches, since a z=2.9676 Lyman limit system cuts out shorter wavelength UV emission, so DLAs cannot be directly 
detected below z=1.976. However, several complex metal line systems have been detected at 0.66<z<2.9 (Giallongo 
et al. 1993, D'odorico et al. 1998). In particular, the system at z=0.6631 has strong low ionization CaII3934-3969 and 
MgII2796-2803 absorption. The Mgll column is 4(±1) x 10^^ cm"^ (D'odorico et al. 1998), which would correspond 
to a hydrogen column of ~ 10^^ cm~^ for solar abundances, assuming that the Mgll ion is the most populated of the 
Mg ions. Since other Mg ions will likely to be populated the hydrogen column could actually be higher than this value. 
Therefore a DLA could well be associated to the 0.6631 system and could be responsible for the X-ray absorption. 
Strong, T > 1, low ionization oxygen edges at ^ 0.33 keV are predicted if N^^ 10^^ cm~^ and if the gas ionization is 
as low, as suggested by the strong Mgll line. Furthermore, at z=2.63-2.82 there are 5 metal systems within a velocity 
range of 800 km/s (Giallongo et al. 1993). They clearly indicate a rich environment at this redshift in the direction of 
PKS2126-158, which could again well be associated with the site of the X-ray absorption. 

On the other hand, a highly ionized absorber at the quasar redshift provides a good fit only if the iron abundance 
is smaller than ~ 0.3 solar, while those of all other elements are fixed to the solar value. The ionized absorber could be 
the resuh of the interaction of the radio jet of PKS2126-158 with the ISM of the host galaxy. PKS2126-158 is a GPS 
quasar and GPS sources are often explained in term of a young, newly born radio sources. Silk & Rees (1998) and 
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Fabian (1999) proposed a scenario in which newly born, high luminosity quasars, are still enshrouded by large densities 
of soft X-ray absorbing gas. If its X-ray absorption is intrinsic PKS2126-158, could then be an example of these newly 
born, highly obscured quasars. The connection of X-ray cut-offs and newly born quasars is highly suggestive, and it 
must be studied in detail with future high resolution X-ray observation, which can be able to constrain tightly the 
redshift of the absorbers. 

In the ionized absorber solution the limit to the iron abundance is rather tight (Z^g < 0.3). This is due to the 
good sensitivity of the LEGS and MEGS at 1.5-2 keV, where the Fe-K edges are redshifted. This limit contrasts with 
the estimate of the nuclear iron abundance obtained by Thompson et al. (1999) using the intensity of iron optical 
emission lines. However, as mentioned in the introduction optical estimates may be inaccurate. The X-ray metallicity 
determination and optical lines measurements could be combined to constrain the line emitting gas physical status. 
An underabundance of iron ([a/Fe]> 0.5), would be consistent with a scenario in which the ISM has been enriched 
mostly by type II SNe. If confirmed, this could be used to put constraints on the time needed to form the bulk of type 
la SNe, i.e the time to form a WD in a binary system plus the time needed to accrete matter from the companion 
to surpass the Ghandrasekar limit, which brings to the SN explosion. All these time-scales are today rather poorly 
defined. 

Since iron aggregates easily in dust grains, an iron abundance of the X-ray absorbing matter significantly less than 
the solar value can contribute to explaining the lack of strong ultraviolet reddening in this source (Elvis et al. 1994, 
Bechtold et al. 1994), expected from the measured X-ray absorbing column density, if the dust to gas ratio and dust 
composition are similar to the Galactic one. A low iron abundance absorber can also help in solving the apparent 
paradox of hard X-ray selected high redshift quasars which appear to be obscured in X-rays but show strong broad 
optical and UV lines and little extinction in the optical band (see Fiore et al. 2001 and Gomastri et al. 2001, Maiolino 
et al. 2001a). Other possible solutions of these contradictions are: (1) the X-ray absorber could be highly ionized 
and the dust destroyed, if the matter is close enough to the active nucleus; (2) the dust grains could have different 
properties with respect to the "standard" Galactic mixture (Maiolino et al. 2001b). 

Although the ionized absorber hypothesis is statistically tenable and astrophysically intriguing, in that it helps in 
resolving the contradiction between the X-ray absorption and the lack of UV extinction, the present low resolution 
data are inadequate to choose this as a unique solution. The precise determination of the absorber redshift, as well as 
of its metallicity and ionization state must await for the high resolution observations that have become possible with 
the gratings on board of the Ghandra and XMM-Newton satellites. 
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